Environmental noise reduction for holonomic quantum gates 



Daniele Parodi, 1,2 Maura Sassetti, 1 ' 3 Paolo Solinas, 4 and Nino Zanghi 1 ' 2 

1 Dipartimento di Fisica, 
Universita di Genova, Genova, Italy 
2 Istituto Nazionale di Fisica Nucleare (Sezione di Genova), 
Genova, Italy 3 INFM-CNR Lamia 
Via Dodecaneso 33, 16146 Genova, Italy 
4 Laboratoire de Physique Theorique de la Matibre Condensee, 
Universite Pierre et Marie Curie, 
Place Jussieu, 75252 Paris Cedex 05, France 
(Dated: February 1, 2008) 

We study the performance of holonomic quantum gates, driven by lasers, under the effect of 
a dissipative environment modeled as a thermal bath of oscillators. We show how to enhance the 
performance of the gates by suitable choice of the loop in the manifold of the controllable parameters 
of the laser. For a simplified, albeit realistic model, we find the surprising result that for a long time 
evolution the performance of the gate (properly estimated in terms of average fidelity) increases. On 
the basis of this result, we compare holonomic gates with the so-called Stimulated Raman adiabatic 
passage (STIRAP) gates. 

PACS numbers: 03.67.Lx 



I. INTRODUCTION 

The major challenge for quantum computation is posed 
by the fact that generically quantum states are very del- 
icate objects quite difficult to control with the required 
accuracy — typically, by means of external driving fields, 
e.g., a laser. The interaction with the many degrees of 
freedom of the environment causes decoherence; more- 
over, errors in processing the information may lead to a 
wrong output state. 

Among the approaches aiming at overcoming these dif- 
ficulties are those for which the quantum gate depends 
very weakly on the details of the dynamics, in particu- 
lar, the holonomic quantum computation (HQC) [l| and 
the so-called Stimulated Raman adiabatic passage (STI- 
RAP) 0,0, In the latter, the gate operator is obtained 
acting on the phase difference of the driving lasers dur- 
ing the evolution, while in the former the same goal is 
achieved by exploiting the non-commutative analogue of 
the Berry phase collected by a quantum state during a 
cyclic evolution. Concrete proposals have been put for- 
ward, for both Abelian [5|, |6[ and non-Abelian holonomies 
0, S, i, E3, El El- The main advantage of the HQC 
is the robustness against noise derivin g fr om a imperfect 
control of the driving fields 0, El, E1E1, E3, El El, HH • 

In a recent paper plj we have shown that the dis- 
turbance of the environment on holonomic gates can be 
suppressed and the performance of the gate optimized 
for particular environments (purely superohmic thermal 
bath). In the present paper we consider a different sort 
of optimization, which is independent of the particular 
nature of the environment. 

By exploiting the full geometrical structure of HQC, 
we show how the performance of a holonomic gate can 
be enhanced by a suitable choice of the loop in the man- 
ifold of the parameters of the external driving field: by 



choosing the optimal loop which minimizes the "error" 
(properly estimated in terms of average fidelity loss). Our 
result is based on the observation that there are different 
loops in the parameter manifold producing the same gate 
and, since decoherence and dissipation crucially depend 
on the dynamics, it is possible to drive the system over 
trajectories which are less perturbed by the noise. For a 
simplified, albeit realistic model, we find the surprising 
result that the error decreases linearly as the gating time 
increases. Thus the disturbance of the environment can 
be drastically reduced. On the basis of this result, we 
compare holonomic gates with the STIRAP gates. 

In Sec. II the model is introduced and the explicit 
expression of the error is derived. In Sec. Ill we find 
the optimal loop, calculate the error, make a comparison 
with other approaches, and briefly sketch how to treat a 
different coupling with the environment. 

II. MODEL 

The physical model is given by three degenerate (or 
quasidegenerate) states, |+), |— ), and |0), optically con- 
nected to another state \G). The system is driven by 
lasers with different frequencies and polarizations, acting 
selectively on the degenerate states. This model describes 
various quantum systems interacting with a laser radia- 
tion, ranging from semiconductor quantum dots, such as 
excitons |l2| and spin-degenerate electron states 0| , to 
trapped ions 0| or neutral atoms 0|. 

The (approximate) Hamiltonian modeling the effect of 
the laser on the system is (for simplicity, h = 1) 0, E2| 

Ho(t)= E m{3\ + (e- ie %(t)\j){G\+H.c)], (1) 

J=+,-,0 

where £lj(t) are the timedependent Rabi frequencies de- 
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pending on controllable parameters, such as the phase 
and intensity of the lasers, and e is the energy of the 
degenerate electron states. The Rabi frequencies are 
modulated within the adiabatic time t ac i, (which coin- 
cides with the gating time), to produce a loop in the pa- 
rameter space and thereby realize the periodic condition 
H (t ad ) = H (Q). 

The Hamiltonian (TT]) has four time dependent eigen- 
states: two eigenstates \Ei{t)) , i = 1,2, called bright 
states, and two eigenstates \Ei(t)), i = 3,4, called dark 
states. The two dark states have degenerate eigen- 
value e and the two bright states have timedependent 
energies \±(t) = [e ± y/e 2 + 4tt 2 (t)}/2 with Q 2 (t) = 

Ei =± ,ol^)l 2 [H- 

The evolution of the state is generated by 



U t = Te~ l ti dt ' Ho{t '\ 



(2) 



where T is the time-ordered operator. In the adiabatic 
approximation, the evolution of the state takes place in 
the degenerate subspace generated by |+), |— }, and |0). 
This approximation allows to separate the dynamic con- 
tribution and the geometric contribution from the evolu- 
tion operator. Expanding Ut in the basis of instantaneous 
eigenstates of H {t) (the bright and dark states), in the 
adiabatic approximation, we have 

u t ^Y, e ~ ip ° Ei{t ' )dt '\ E ^ E ^)\ w *> ( 3 ) 



where 



U t =Te^ dTV ^\ (4) 



here V is the operator with matrix elements Vij (t) — 
(Ei(t)\dt\Ej(t)). The unitary operator Ut plays the role 
of timedependent holonomic operator and is the funda- 
mental ingredient for realizing complex geometric trans- 
formation whereas X) - e _i /o E J^ d ''|£|,-(t)}{JSj(t)| is the 
dynamic contribution. 

Consider U t for a closed loop, i.e., for t — t a 4, 



U = U, 



(5) 



If the initial state |^o) is a superposition of |+) and |— ), 
then U\if>o) is still a superposition of the same vectors (in 
general, with different coefficients) [l2T] . Thus the space 
spanned by |+) and — ) can be regarded as the "logi- 
cal space" on which the "logical operator" U acts as a 
"quantum gate" operator. Note that for t < t a d, Ut\ipo) 
has, in general, also a component along |0). However, 
as it is easy to show [22|, at any instant t < t a a, Ut\4>o) 
can be expanded in the twodimensional space spanned 
by the dark states \Es(t)) and \Ei(t)). It is important 
to observe that U depends only on global geometric fea- 
tures of the path in the parameter manifold and not on 
the details of the dynamical evolution [l|, [l2j • 

To construct a complete set of holonomic quantum 
gates, it is sufficient to restrict the Rabi frequencies 



£lj(t) in such a way that the norm of the vector 
£1 = [Clo(t), Q_ (£)] is time independent and the 

vector lies on a real three dimensional sphere [1, H3 |. 
We parametrize the evolution on this sphere as f2+(t) = 
sin 6>(i)cos (f>(t), = sin 0(*)sin (f>(t) and fl (t) = 

cos 8(t) with fixed initial (and final) point in 6(0) = 0, 
the north pole By straightforward calculation we obtain 
the analytical expression for V(t) in eq. (fj|, V(t) — 
icr y cos[9(t)]<j)(t), where a y is the usual Pauli matrix 
written in the basis of dark states. Thus, the oper- 
ator becomes Ut — cos[a(t)] — ia y sin[a(i)], here 
a(t) = J Q dr0(r) cos 9(t). Accordingly, the logical op- 
erator U ([5]) is 



where 



U 



a(t ad ) 



ut,, sm a, 



dr^(r) cos 9(t) 



(6) 



(7) 



is the solid angle spanned on the sphere during the evolu- 
tion. Note that the are many paths on the sphere which 
generate the same logical operator U, and span the same 
solid angle a. 

In a previous work we have studied how interaction 
with the environment disturbs the logical operator U [2l| . 
The goal of the present paper is to analyze whether and 
how such a disturbance can be minimized for a given 
U. To this end, we model the environment as a thermal 
bath of harmonic oscillators with linear coupling between 
system and environment [23j. The total Hamiltonian is 
then 



iY 



H = H a (t) + J2( 



2m n 



^m a ulx 2 a 



+ c a x a A), (8) 



where A is the system interaction operator called, from 
now on, noise operator. 

We now consider the time evolution of the reduced 
density matrix of the system, determined by the Hamil- 
tonian ©. We rely on the standard methods of the "mas- 
ter equation approach," with the environment treated in 
the Born approximation and assumed to be at each time 
in its own thermal equilibrium state at temperature T. 
This allows to include the effect of the environment in 
the correlation function (ks = 1) 



g(r) = / J(«) 



roth ( — — j cos(o;t) —i sin(wr) 



Here the spectral density is 



JV 



2 ' m a uj a 

a—l 



8(u - LU a ), 



du. 

(9) 

(10) 



at the low frequencies regimes, is proportional to io s , with 
s > 0, i.e., s = 1 describes a Ohmic environments, typ- 
ical of baths of conduction electrons, s = 3 describes 
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a super- Ohmic environment, typical of baths of phonons 
[2l|,[24j]. The asymptotic decay of the real part of g(r) de- 
fines the characteristic memory time of the environment. 
Denoting with p(t) the time evolution of the reduced den- 
sity matrix of the system in the interaction picture, e.g., 
p(t) = ujpUt, one has [H] 



their degeneracy, shifting one of them. In spite of its sim- 
ple form, this A is nevertheless a realistic noise operator 
for physical semiconductor systems 0]. 



III. MINIMIZING THE ERROR 



fttad) = P(0) + 



dt / dr{g(T)[AA' p(t - r) - A'p(t - r)A] 



(11) 



+ g(-r)[p(t~r)A'A-Ap(t-r)A']. 



Here A and A' stand for A(t) and A(t — r), with the tilde 
denoting the time evolution in the interaction picture. 

In quantum information the quality of a gate is usually 
evaluated by the fidelity T, which measures the closeness 
between the unperturbed state and the final state, 



T= Gfo(0)|wV(i o< ,)W|M0)), 



(12) 



where |^o(0)} is the initial state, and p(t a <i) = Up(t a d)W 
is the reduced density matrix in the Schrodinger picture 
starting from the initial condition p(Q) = \ipo(0)) (ipo(0)\. 
The average error is defined as the average fidelity loss, 
i.e., 

s =< 1 - F>= l- < (Mo)\p(t ad )\Mo)) >, (13) 

where < • • • > denotes averaging with respect to the 
uniform distribution over the initial state \ipo(0)). 

The right-handside of Eq. (fl~3"|) can be computed by 
the following steps: 

(1) solving Eq. pip in strictly second order approxima- 
tion; this approximation corresponds to replace p(t — r) 
with p(0); 

(2) using the adiabatic approximation U(t — r,t) « 
exp(rri?o(t)); 

(3) expanding the scalar product in Eq. (|13|) with respect 
to a complete orthonormal basis {\<p n (t))} 7 n — 1,2,3, 
orthogonal to \ipo(t)}. In this way, one obtains 

S = ^J* ad dt G(t)\{Mt)\A\Mt))\ 2 ^ (14) 
where 



G(t) = / dr{Re[5(r)] cos(wont) + Im[g(r)] sin(cjont))}- 
Jo 

(15) 

Here, ajQ n = luq — u> n are the energy differences associated 
to the transition ipo «-> </>„, with luq = e, lj\ = A+, loi — 
A_, and W3 = e. 

The interaction between system and environment is ex- 
pressed by the noise operator A in Eq. ||5J). We shall now 
make the assumption that A = diag{0, 0, 0, 1} in the \G), 
±), and |0) basis. In this case the transition between de- 
generate states are forbidden, however the noise breaks 



The problem can be stated in the following way: given 
the noise operator A and the logical operator U, find a 
path on the parameter space (the surface of the sphere, 
described above) which minimizes the error 5. 

The total error S, given by Eq. (TT4")) , can be decom- 
posed as 



S = 5 tr + S 



pd • 



(16) 



where the transition error, 5 t r, is the contribution to the 
sum of the nondegenerate states (coo n 7^ 0) and the pure 
dephasing error S p d is the contribution of the degenerate 
states (u)Q n — 0). Thus 



x n 

Ond — — 



dt 



, J(oj) . I u> 

aid coth — 

2T 



1 



and 
8tr : 



E 



sin(wf) (1 + 1 sin 2 2a(t) ) sin 4 6(t) (17) 



r „ / sin 2 26{t)dt, 


(18) 



where 

Ton = >/(|wOn|) 



coth 



w 0n 



2T 



- sgn(won) 



(19) 



correspond to the transition rates calculated by standard 
Fermi golden rules, supposing, as usual, G(t) ~ G(oo) for 
g(j) strongly peaked around r = 0. In the following we 
define for simplicity 



1 



On- 



Since we are interested at long time evolution, we start 
discussing the transition error which dominates in this 
regime [111]. 



A. Transition rate 

As explained in Sec. [Til the holonomic paths are closed 
curves on the surface of the sphere which start from the 
north pole. It turns out that the curve minimizing 5t r 
can be found among the loops which are composed by a 
simple sequence of three paths (see the Appendix): evo- 
lution along a meridian (0 = const), evolution along a 
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FIG. 1: The error Str versus 8m for two different a values: 
a — 7r/2 (dashed line) and a — n/4 (full line) correspond to 
NOT and Hadamard gate, respectively. 

parallel {9 = const) and a final evolution along a merid- 
ian to come back to the north pole. 

The error S tr in (fTS]) , depends on a given by Eq. |J7J, 
9m (the maximum angle spanned during the evolution 
along the meridian), A<fi (the angle spanned along the 
parallel), and angular velocity v. We allow A<f) > 2tt 
which corresponds to cover more than one loop along 
the parallel. The velocity along the parallel is v(t) — 
4>{t) sin 9 and that along the meridian is v(t) = 9{t). In 
the following we assume that v is constant, and it cannot 
exceed the maximal value of i>max, fixed by adiabatic 
condition i>max <C fl- 

The parameters a, 9m, and A<ji> are connected by the 
relation a — A</>(1 — cos 9m)- The error 5t r is then 

8 tr = 6% + 6? r , (20) 

where 

5g = *(e M -\ftoi4Bj\ (21) 

is the contribution along the meridian and 
P r ,asin 9 M sin 2 29 M 

°tr — K ; n 

v 1 - cos 9 M 

is the contribution along the parallel. 

In Fig.[T]<5i r is plotted for a — it/2 and a = 7r/4 (corre- 
sponding to NOT and Hadamard gate, respectively) as a 
function of 9m- One can see that 5t r has a local minimum 
for 9m = 7r/2 and a global minimum for 9m — where 
the error vanishes. This suggests that the best choice is 
to take 9m as small as possible. 

It is interesting to consider the dependence of St r also 
on the evolution time t a d- For simplicity, we set the ve- 
locity v = Umax- In this case, changing 9m (and then 
A(f>) corresponds to a change in the evolution time. We 
obtain 

9 M = arccos ( 1 - - J , (23) 



FIG. 2: The error Str versus umaxiad for two different a val- 
ues: a = 7r/2 (dashed line) and a = n/4 (full line) correspond 
to NOT and Hadamard gate, respectively. The dotted-dashed 
line shows the value of the error at 6 = 7r/2. The circles show 
the critical value of wmaxtod above which the best loop is the 
one with the minimal 6m- 



where 

m = i^T [( Umaxtad ) 2 + ° 2 ] ■ ( 24 ) 

Using these relations, 5^ and 5f r , given by (f2"Tj) and 
([22)1 become functions of t a d, fmax, and a. Note that 
m measures the space covered along the parallel, in fact 
A</> = 27rm. 

In Fig. [5] we see the behavior of S tr as a function 
of fmax^ad- The first minimum for both curves corre- 
sponds to 9m = ?r/2, then the curves for long t a d de- 
crease asymptotically to zero corresponding to the region 
in which 9m — * 0. In this regime we have St r oc l/t a d 
which is drastically different from the results obtained 
with other methods where 5t r oc t a d, (see Refs [U, |2o| 
and below Sec. IIII Cp . It should be observed that this 
surprising results is a merit of holonomic approach which 
allows to choose the loop in the parameter space, with- 
out changing the logical operation as long as it subtends 
the same solid angle. Observe that small 9m and long 
t a d mean large value of to, i.e., multiple loops around the 
north pole. 

Figure [3 shows that, for a given gate, there is a criti- 
cal value k c of Vmaxiad which discriminate between the 
choice of 9m (e.g., k = 6 for the Hadamard gate and 
k = 25 for the NOT gate). For vm&xtad < k c the best 
choice for the loop is 9m = tt/2; For fmax^ad > k c the 
best choice is the value of 9m determined by eq. ([23]) and 

(2D. 

Note that the region Vm&xt a d > k c is accessible with 
physical realistic parameters [T^ |. For example, if wc 
choose the laser intensity SI = 20 meV and Wmax = S^/50 
(for which values the nonadiabatic transitions are forbid- 
den), the critical parameter corresponds to the critical 
time of 15 ps for the Hadamard gate and 42 ps for the 
NOT gate. 
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B. Pure Dephasing 

Until now we have ignored the pure dephasing effect 
because we have assumed that it is negligible in com- 
parison with the transition error for long evolution time. 
Now, we check that the pure dephasing error contribu- 
tion can indeed be neglected. We can write the pure 
dephasing error using Eq. (| 1 T() and splitting to parallel 
and meridian part as 



°pd 



dt 



au — 1 — coth 



-) 

2TJ 



Q[a(t)] sin cot sin 9 



M 



(25) 



and 



M _ i "max 

°pd — 



dt 



JO 

4/ 



sin uiti sin (fmaxi) + sin 



1 - 



where Q[a(t)] = 1 + 1/2 sin 2 [2a(t)]. 

To estimate S p d we assume that t a d is longer with re- 
spect to the characteristic time of the bath. Remember- 
ing that J(co) oc lu s , the pure dephasing error behavior 
along the parallel part at the temperature T is 



to the north pole along meridians. Since the loop, as 
in our model, is a sequence of meridian-parallel-meridian 
path, we can calculate the error and make a direct com- 
parison. In this case, the transition error results propor- 
tional to 6tr oc t ac i and grows linearly in time while for 
HQC Str oc 1 /t a d- Therefore, the HQC is fundamentally 
the favorite for long application times with respect to the 
STIRAP ones. 

Moreover, we can show that the freedom in the choice 
of the loop allows us to construct HQC which perform 
better than the best STIRAP gates. In Ref. Q the mini- 
mum error (not depending on the evolution time) for STI- 
RAP was obtained reaching a compromise between the 
necessity to minimize the transition, pure dephasing error 
and the constraint of adiabatic evolution. With realistic 
physical parameters [2l[ (J(ui) = koj^e^/"") 2 ; n = 10 



meV, e = leV, v„ 



fi/50, k = lQ- 2 {meV)- 2 , lu c = 0.5 



meV and for low temperature), the total minimum error 

With the same parame- 



= 10 



-3 



in Ref. [4j is <5 stirap 
ters, we still have the possibility to increase the evolution 
time in order to reduce the environmental error. How- 
ever, for evolution time t a d = 50 ps we obtain a total 
error 5 = 1.5 X 10~ 4 for the NOT gate and 8 = 4 x 10~ 5 
for the Hadamard gate, respectively. As can be seen, the 
logical gate performance is greatly increased. 



D. More general noise 



s+3 



K fLs+2 (27) 
,T»l/tad 



while the along meridian is 



tad 



(28) 



Then, we can conclude that the pure dephasing can al- 
ways be neglected for long time evolution because it de- 
creases faster than the transition error. 



C. Comparison between HQC and STIRAP 

We make a comparison between holonomic quantum 
computation (HQC) and the STIRAP procedure which is 
an analogous approach to process quantum information. 
The STIRAP procedure ([2|, |4|) is, in its basic points, 
very similar to the holonomic information manipulation. 
The level spectrum, the information encoding, the evolu- 
tion produced by adiabatic evolving laser are exactly the 
same. The fundamental difference is that in STIRAP 
the dynamical evolution is fixed (we must pass through 
a precise sequence of states) and then the correspond- 
ing loop in the parameter space is fixed. In particular, 
we go from the north pole to the south pole and back 



Until now we have discussed the possibility to minimize 
the environmental error by choosing a particular loop in 
the parameter sphere but the structure of the error func- 
tional clearly depends on the system-environment inter- 
action. Then one might wonder if the same approach can 
be used for a different noise environment. 

For this reason, we now briefly analyze the case of noise 
matrix in the form A = diag{0, 1, 0, —1}. Again, for long 
evolution we can neglect the contribution of the pure de- 
phasing and focus on the transition error. In this case the 
interesting part of the error functional takes the form 



6 tr = if [(^ sin 26» cos 26) 2 + (sin 9 sin 2<j)) 2 ]. (29) 

Even if the analysis in this case is much more com- 
plicated, it can be seen that 5t r has an absolute mini- 
mum for 9m = 0. The long time behavior is the same 
(Str oc 1 /t a d) such that the results are qualitatively anal- 
ogous to the above ones: for small 9m loops (or long 
evolution at fixed velocity) the holonomic quantum gate 
presents a decreasing error. Then even in this case it is 
possible to minimize the environmental error. 



IV. CONCLUSIONS 

In summary, we have analyzed the performance of 
holonomic quantum gates in the presence of environmen- 
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tal noise by focusing on the possibility to have small 
errors choosing different loops in the parameter mani- 
fold. Due to the geometric dependence, we can imple- 
ment the same logical gate with different loops. Since 
different loops correspond to different dynamical evolu- 
tions, we have used this freedom to construct an evolu- 
tion through "protected" or "weakly influenced" states 
leading to good holonomic quantum gates performances. 
This allows to select (once that the physical parameter 
are fixed) the best loop which minimizes the environ- 
mental effect. (Note that this optimization procedure 
is rather independent of the details of the simple model 
we have considered and arguably, it could be extended 
to more complicated systems without any substantial 
modification.) We have shown that for long time evo- 
lutions the noise decreases as l/t a d while in the other 
cases it increases linearly with adiabatic time. We also 
have shown that the same features can be found with 
different kinds of noise suggesting the possibility to find 
a way to minimize the environmental effect in the pres- 
ence of any noise. These results open a new possibility 
for implementation of holonomic quantum gates to build 
quantum computation because they seem robust against 
both control error and environmental noise. 
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to C n can approximate any closed curve on the sphere, 
the best curve can be found in C\. 

The crucial point is to show that any curve in C2 can be 
replaced by a curve in C\ . Let us consider a generic curve 
C2 in C2 spanning a solid angle a: composed by a seg- 
ment of a meridian (with 9 going from to 9i), a parallel 
(spanning a A0i angle), meridian (with 9 : 9\ — > 62), a 
parallel (spanning a A</> 2 angle), and finally a segment 
to the north pole along a meridian. Let us consider two 
closed curves C\ and C\ in C\ subtending the same solid 
angle a with, respectively, 9\ and 9i as maximum angle 
spanned during the evolution along the meridian. First 
we analyze (j2"0)l along the meridian. Without losing gen- 
erality, we can take 9\ < 6*2 ; it is clear from Eq. ([2"Tj) that 
the value of 5t r along the meridian for C\ is smaller that 
for C\: < S^{. We note from the Eq. (T2T]), suitable 

extended to C 2 , that the two paths along the meridians 
depends only on #2 and then produce the same error of 

= (Al) 

The difference between the contribution along the par- 
allel is 



6Z -SZ 2 = A0i ( sin 6»i sin 2 29 x sin 9 2 sin 2 29 2 

1 \ 1 - cos 9 2 

(A2) 

and 



APPENDIX A: MINIMIZING THEOREM 

Let us consider the family C n composed of the closed 
curves generated by a sequence of n paths along a parallel 
(9 = const) alternated with paths along a meridian (0 — 
const). We call C n a generic curve in this family. For 
example, the family C\ contains all the closed curves com- 
posed by the sequence of path meridian-parallel-meridian 
while the family C 2 contains the curves meridian-parallel- 
meridian-parallcl- meridian. 

We argue that the closed curve minimizing the error 
in Eq. (fT5|) can be found in the C\ family. First, we 
show that any closed curve in C2 spanning a solid angle 
a on the sphere can be replaced by a closed curve in C\ 
spanning the same angle and producing a smaller error. 
In analogous way any closed curve in C3 can be replaced 
by a closed curve in C2 with smaller error and so on. By 
induction we obtain that any closed curve in C n can be 
replaced by a curve in C\ spanning the same solid angle 
but producing smaller error. Since the curve belonging 



&p -8^ = Afoi sin 6» 2 sin 2 29 2 # sin 9 X sin 2 26»i 

2 1 \ 1 - COS 9\ 

(A3) 

Analysis of the positivity of the quantities given by 
Eqs. (|A2|) and (|A3|) shows that S^ 2 cannot be at the 
same time smaller than 6£ 2 and 8^ 2 . In fact, there are 
two possibilities: If S£ 2 > 5^, from Eq. (jATj) and (K3\t , 

5c 2 = S% + C > 5c\ + 5 c } = 5 c{ , (A4) 

and the best closed curve is C\. If 8^ 2 > S^ 2 , from Eqs. 
(IAT1) and fA2]l . 

Sc a = S&+Sg a >6% + S* f = S c , i (A5) 

and the best closed curve is C 2 . 

In the same way it can be shown that any closed curve 
in C3 can be replaced by a closed curve in C2 with smaller 
error. 
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